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ABSTRACT 

We analyze the possible existence of a pulsational instability excited by the 
e-mechanism during the last few centuries of evolution of a Chandrasekhar mass 
white dwarf prior to its explosion as a Type la supernova. Our analysis is moti- 
vated by the temperature sensitivity of the nuclear energy generation rate (~ T 23 ) 
in a white dwarf whose structural adiabatic index is near 4/3. Based upon a lin- 
ear stability analysis, we find that the fundamental mode and higher order radial 
modes are indeed unstable and that the fundamental mode has the shortest 
growth time scale. However, the growth time scale for such instability never be- 
comes shorter than the evolutionary timescale. Therefore, even though the star 
is pulsationally unstable, we do not expect these radial modes to have time to 
grow and to affect the structure and explosion properties of Type la supernovae. 



Subject headings: stars: supernovae — stars: pulsations 



-2- 



1. Introduction 

The explosion of a Chandrasekhar mass white dwarf as a Type la supernova is one of the 
most dramatic events that can befall a star and, at the same time, one of the most complex. 
The outcome depends on details of the ignition process (Woosley, Wunsch, & Kuhlen 2004) 
and the uncertain physics of turbulent flame propagation (see Hillebrandt & Niemeyer 2000, 
Gamezo et al. 2003, and references therein). Both of these depend on the presupernova 
evolution of the star which sets e.g., the central density at the time of runaway and the 
distribution of ignition points. However, since pioneering studies in the 1970's by Paczynski 
(1972), Couch & Arnett (1973, 1974, 1975), and Iben (1978), comparatively little attention 
has been given to the evolution of the interior of the presupernova star. It is known that 
the star ignites carbon centuries before it finally explodes, commencing a long ramp up that 
eventually concludes when convection can no longer carry away the excess energy generated 
by the degenerate nuclear runaway. It is during this phase that the URCA process may 
lead to a complex convective structure in the core and perhaps even oscillations related to 
the switching on and off of URCA neutrino losses (Iben 1978). In this paper we consider a 
comparatively simple question: "Is presupernova carbon burning inherently unstable - even 
in the absence of URCA losses" ? 

Carbon burning occurs at a rate that is extremely temperature sensitive, at least T 23 
even at the end of the stable burning near T « 7x 10 8 K. The white dwarf is also supported 
chiefly by degenerate electrons that are relativistic, hence T « 4/3. In other contexts in 
stellar evolution, this combination often tends towards instability (e.g., Ledoux 1941; Maeder 
1985; Baraffe, Heger, & Woosley 2001). Large excursions in the radius are not strongly 
damped and the non-linearity of the reaction rate promotes instability. This is known as 
the e-mechanism (Unno et al. 1989). Such an instability leads to oscillations around the 
hydrostatic equilibrium configuration and periodic phases of expansion and contraction on 
a dynamical timescale Td yn ~ (Gp)^ 1 ^ 2 . These types of oscillations are described as acoustic 
modes. They can be excited if an excitation mechanism operates in regions where the 
amplitude of eigenfunctions is large. 

We have thus undertaken a linear stability analysis of the structure of Chandrasekhar 
mass white dwarf models prior to explosion, based on a similar analysis performed for very 
massive stars (Baraffe et al. 2001). We pick up the white dwarf when it has just ignited 
carbon (nuclear energy generation substantially in excess of plasma neutrino losses) and is 
developing an extended convective core. The central temperature at this point is ~ 3 x 10 8 K 
and it is about 300 years until the star's death. The evolutionary models and input physics 
used are described in §2 and the results of the linear stability analysis are presented in §3. 
A discussion follows in 84. 
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2. Evolutionary models prior to explosion 



2.1. Characteristic evolutionary timescales 

Following ignition, a convective region develops in the white dwarf that eventually grows 
to encompass most of its mass. Heat is transported outwards, but both radiative losses and 
neutrino losses are negligible. Instead, the energy mostly goes into heating the convective 
portion of the star, and, to a lesser extent, expansion. The heat capacity is given by 
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(Woosley et al. 2004), where e is the specific internal energy, T 8 = T/10 8 K and pg = 
p/10 9 gcm -3 . For the relevant temperature and density range, T 8 = 2 to 7, pg = 1 to 3, 
the ionic term dominates with a minor contribution from the electrons. To a factor-of-two 
accuracy, cp ~ 10 15 erg g" 1 (10 8 K) _1 . The total thermal energy is 
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The correct definition for pure "thermal energy" should be expressed in terms of cy 
rather than cp, since cp also contains an expansion term PdV . However, because of strong 
degeneracy in the interior of the white dwarf, cy ~ cp and Eq. (2) provides a good estimate 
of the thermal heat content. Moreover, for our purpose below, to estimate the gravothermal 
heat content, this is the appropriate quantity to use. If the mass-averaged temperature is 
about half the central value and the convection zone, about one solar mass, 
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where the index "c" stands for central quantities of the star. Actually, from computer models 
discussed in the next section we know that the convective core grows from 0.2 M Q to 1.15 
M Q as the central temperature rises from 3 x 10 8 K to 7 x 10 8 K (see Table 1). A better 
estimate of the heat in the convective region is 
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This actually agrees to much better than factor-of-two accuracy with the detailed computer 
results described in the next section. 
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The total nuclear power generated by the star is (Woosley et al. 2004) 




(5) 



where v ~ 23 for the relevant range of temperature and X\i is the carbon mass fraction. 
This expression is valid in the temperature range 4 < T 8 < 7 and becomes increasingly 
inaccurate at lower temperature because of temperature dependence of v and especially the 
complex dependence of the screening correction on temperature and density. 

The nuclear timescale, which characterizes the increase of nuclear energy generation 
rate in the center, taking into account the reaction of the rest of the star, is 



where it is assumed that L = dH/dt (see Table 1). 

Since the central temperature, or the total heat content of the star need only increase 
~3% for the power to double, the white dwarf evolution is governed by the rise of the power, 
L. Thus, the nuclear timescale tl and the time to explosion are the same and both are 



The star takes a few centuries to go from ignition at 3.0 x 10 8 K to 7 x 10 8 K, but explodes a 
few minutes later. This is a considerably longer time scale than one would get by evaluating 
a local nuclear timescale r nuc = (cpT) / (ue nuc ) , just from the central energy generation rate 
and heat capacity, because the cooler middle regions of the white dwarf act as an appreciable 
heat reservoir for the energy generation. Until the central temperature exceeds 7 x 10 s K, 
convection keeps the two efficiently coupled (see discussion in §4). 
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considerably shorter, by about a factor of 5, than the thermal time scale th 
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Table 1: Nuclear and explosion timescales (see §2) 



T 8 , c 


P9,c 


M-conv 
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dH 
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^expl 






(M©) 


(erg/s) 


(erg/s) 


(y r ) 


(y r ) 


3.0 


3.60 


0.2 


2.7xl0 36 


2.3xl0 36 


210 


300 


4.0 


3.50 


0.6 


2.1xl0 39 


1.9xl0 39 


0.75 


1.5 


5.0 


3.35 


0.85 


3.3xl0 41 


2.9xl0 41 


0.010 


0.017 


6.0 


3.16 


1.05 


1.9xl0 43 


1.6xl0 43 


3.5xl0~ 4 


4.6 xl0~ 4 



NOTE: T 8 c is the central temperature in units of 10 s K; p 9 c is the central density in units 
of 10 gem ; M conv the convective core; L the nuclear power defined in Eq. (5); ^ the 
time derivative of the total heat content H defined in Eq. (4); r L and t exp i are respectively 
the nuclear and the explosion timescales in years. 
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2.2. Stellar Models 

Computer models of accreting white dwarfs approaching criticality were constructed 
using the KEPLER stellar evolution code (Weaver et al. 1978). A composition of 30% 
12 C and 70% 16 0, by mass fraction was assumed, though the outcome will not depend 
appreciably on this assumption. The initial model consisted of a "warm" white dwarf with 
central temperature 10 8 K and a mass of 2.6 x 10 33 g (~ 1.307 M Q ). Test calculations that 
started from a cooler white dwarf resulted in a higher mass at ignition, but the results of 
the pulsation analysis were not affected by this modification. The white dwarf then accreted 
matter (also carbon and oxygen in the same proportion) at 10" 7 M Q yr" 1 . The accretion was 
stopped when the runaway had started, but about 80 centuries before the final incineration, 
to allow the addition of the well-resolved surface layers. At this point the star has reached a 
central density of pg jC = 3.64 and a central temperature of Tg ;C = 2.43. In the remaining time 
the star would have accreted less than 1O _3 M , which does not affect the structure of the 
star or the runaway (note that the runaway occurs many times that mass before reaching 
the Chandrasekhar mass). 

The Lagrangian grid uses zone masses ranging from 10 29 g in the center over some 10 31 g 
in the middle of the star to a smooth gradient in zone size down to 10 15 g at the surface. 
Convection is treated using standard mixing length theory with a mixing length L mix = Hp. 

During the last ~ 300 yr prior to explosion the nuclear energy generation rate, e nuc , 
rises rapidly in the center of the white dwarf. This is illustrated in Fig. 1 that displays 
the evolution of central temperature, density and nuclear energy as a function of the time 
remaining before explosion (t exp i = corresponds to the time of flame ignition/ start of 
explosion). While only a small fraction of the energy released goes into expanding the star, 
the expansion is significant because r ~ 4/3. 

3. Pulsation analysis 

The system of equations which are linearized are the basic hydrodynamic equations 
written under the assumption of spherical symmetry (cf. Unno et al. 1989, p. 89). Details 
of the method and pulsation code can be found in Baraffe et al. (2001). The unperturbed, 
equilibrium state is described by hydrostatic equilibrium. The nuclear energy term appearing 
in the energy conservation equation includes neutrino energy loss. An important uncertainty 
in the present analysis is due to the assumption that convection is frozen in, neglecting the 
perturbation of the convective flux. This simplification is based on the argument that the 
convective timescales in the interior of the white dwarf remain significantly larger than the 
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period of pulsation of the excited modes (see below). A sequence of evolutionary models 
is analyzed, starting from a central temperature T 8 c = 3, ~ 300 yr before explosion, to 
^8,c = 7, corresponding to ~ 10 minutes before explosion. All models analyzed show that 
the fundamental and higher order modes are unstable. As in Baraffe et al. (2001), the 
assumed time dependence for the eigenfunctions is of the form exp(iat) x exp(JTt), where 
a = 27r/Il is the eigen-frequency, II the corresponding period and K the stability coefficient. 
The e-folding time characterizing the growth timescale of the pulsation amplitude is defined 
by Td = 1/K. In all models analyzed, this fundamental growth timescale is, by at least one 
order of magnitude, shorter than Td of the higher order modes. Our analysis below can thus 
be restricted to the fundamental mode, since it is the mode which grows the fastest. 

Table 2 summarizes the stability analysis results. The fundamental mode pulsation 
period n is ~ 2 s, close to the dynamical timescale. Note that the shortest convective 
timescale is ~ 10 4 s in the first model analysed (Tg jC = 3) and ~ 100 s in the last models. For 
a model with Tg jC = 5 (~ 6 days before explosion), Fig. 2 displays the work dW/dM, which 
is proportional to the energy transferred to the pulsation, with W the so-called work integral 
(Cox 1980; Unno et al. 1989). Positive work indicates driving regions whereas negative work 
corresponds to damping zones. The positive value of dW/dM in the center is due to the 
perturbation of the nuclear energy generation rate and characterizes the e-mechanism. As 
expected for a gas with Ti ~ 4/3, the displacement 5r/r is essentially constant through the 
whole structure. Damping processes are small in the rest of the white dwarf, and the total 
work integral W is essentially determined by the driving zone in the center. 
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Table 2: Linear stability analysis results for an evolutionary sequence of SN la progenitors 
prior to explosion. 



T 8 , c 


P9,c 


R (Mm) 


^cxpl (s) 


Time \p) 


T d (s) 


n (s) 


3 


3.60 


1.75 


9.68xl0 9 


2.99xl0 8 


2.95xlO n 


1.95 


4 


3.50 


1.77 


4.48xl0 7 


8.31xl0 5 


6.19xl0 8 


1.95 


5 


3.35 


1.79 


5.46xl0 5 


8.63 xlO 3 


4.92xl0 6 


1.96 


6 


3.16 


1.82 


1.47xl0 4 


2.26xl0 2 


1.03xl0 5 


1.97 


7 


2.97 


1.85 


6.57xl0 2 


11.3 


4.22xl0 3 


1.98 



NOTE: Models are characterized by the central temperature T 8 c , central density pg jC and 
white dwarf radius R. t cxp \ is the remaining time before explosion and r nuc = cpT/(ve nnc ) is 
the local nuclear timescale defined in the center. Growth timescale (see text) and period of 
the fundamental mode are denoted respectively by r d and n . 
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4. Discussion 

Linear stability analysis alone does not provide information about the maximum pul- 
sation amplitude which can be finally reached. For this purpose, full hydro dynamical cal- 
culations are needed. However, comparison between the growth timescale of an excited 
mode and the evolutionary timescale t exp \ can indicate whether the perturbation may have 
time to grow in order to reach large amplitudes. Before doing such comparisons, it is useful 
to consider the uncertainties affecting the evolutionary timescale. 

As previously discussed in §2.1, the relevant timescale characterizing the evolution to 
explosion is considerably longer than the local nuclear timescale because of the heat reser- 
voir provided by convective coupling to the rest of the star. For an e-mechanism, the growth 
timescale is expected to be closely related to the local nuclear timescale, since r nuc char- 
acterizes the increase of local nuclear energy generation rate in the central region where 
the e-mechanism takes place (see Fig. 2). The discussion below will thus be based on r nuc 
and its comparison with the time to explosion. Qualitatively, convection can lengthen the 
evolutionary timescale compared to the local nuclear timescale, and thus r nuc < t cxp i ~ 7~l. 

For the sake of a quantitative comparison, we have performed test evolutionary se- 
quences with different convection efficiencies, in terms of the mixing length parameter L mix . 
All test sequences start from the same model with central temperature T 8jC = 3. The nuclear 
timescale r nuc as a function of the remaining time before explosion is displayed in Fig. 3 for 
four different cases of convection efficiency: the standard case with L m ; x = Hp (solid line), 
L m i x = 100-ffp (dashed line), L mix = ifp/100 (dotted line) and a case with no convection 
(dash-dotted line). The long-dashed line corresponds to the case t exp i = T nuc- As expected, 
when convection is inhibited, the runaway timescale is essentially determined by the local 
nuclear timescale. On the opposite, the more efficient the convection, the longer the evolu- 
tionary timescale compared to r nuc . Fig. 3 indicates that even in the unrealistic case with 
-^mix = lOOifp, the evolutionary timescale is at most ~ 70 times greater than r nuc . There- 
fore, although the treatment of convection is crucial and still very uncertain during these 
last phases of evolution, one can reasonably expect that such uncertainty will not lengthen 
the evolutionary timescale by more than 100 x r nuc . 

Fig. 4 displays the growth timescale of the fundamental mode and the local nuclear 
timescale r nuc as a function of the evolutionary timescale t exp \, for the standard case of 
convection efficiency (see also Table 2). As expected and illustrated in Fig. 4, r d follows r nuc , 
with r d /r nuc ~ 5 10 2 - 10 3 , the same relation of proportionality being found independently 
of convection efficiency. Since the tests on convection efficiency suggest that t CX pi < 100 r nuc , 
the quantity thus remains systematically larger than the evolutionary timescale t exp i, as 
illustrated in Fig. 4. For standard convection, t& is larger than t cxp \ by a factor ~ 30 at 
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the beginning of the sequence (Tg jC = 3) down to a minimum factor of ~ 6.5 for Tg >c = 7. 
Similar behavior is found for the test sequences with varying convection efficiency, with the 
smallest ratio r d /t cxp i ~ 5 found for the case L mix = lOOifp. Such results suggest that the 
present vibrational instability does not have time to grow to reach significant amplitudes. 

5. Conclusion 

Based on a linear stability analysis of Chandrasekhar-mass white dwarf models prior 
to explosion, we have found that the central conditions are favorable to the excitation of 
radial eigenmodes through the e-mechanism. Such pulsational instability could in principle 
result in phases of expansion and contraction with growing amplitude and thus affect the 
structure of the supernova progenitor. We also found, however, that the growth timescale 
for the pulsation amplitude remains systematically longer than the evolutionary timescale. 
Such instability does not have time to grow significantly before the explosion of the white 
dwarf. 
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Fig. 1. — Evolution of the central temperature, T c , (in K; solid line), up to 7 10 8 K, central 
density p c (in g cm -3 ; dash-dotted line) and central nuclear energy e nuc (in erg/g/s; dashed 
line) of the white dwarf as a function of the remaining time before explosion t cxp i (in s). 
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Fig. 2. — Differential work dW/dM, in arbitrary units, as a function of mass (in units of 
the total mass) in the interior structure of a model 6 days before explosion (Tg jC = 5). The 
dashed line corresponds to the nuclear energy generation log e nuc (erg/g/s). 
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Fig. 3. — Local nuclear timescale r nuc (see text) as a function of the remaining time before 
explosion t cxp i (in s) for different cases of convection efficiency: L m j x = Hp (solid line), 
-^mix = lOOHp (dashed line, almost indistinguishable from the solid line), L mix = Hp/ 100 
(dotted line) and L mix = (dash-dotted line). The long-dashed line indicates the case where 

7"nuc ^expl- 
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Fig. 4. — Characteristic timescales (in s) as a function of the remaining time before explosion 
t cxp i (in s): growth timescale for the vibrational instability t& (solid line) and local nuclear 
timescale r nuc (short-dashed line). The long-dashed line indicates the ideal case where the 
characteristic timescales are equal to £ exp i. 



